home *** CD-ROM | disk | FTP | other *** search
/ BCI NET / BCI NET Dec 94.iso / archives / grafix / misc / ami2d_2_0.lha / ami2d_dist / ami2d / rexx / init.rexx < prev    next >
Encoding:
OS/2 REXX Batch file  |  1994-09-19  |  2.0 KB  |  119 lines
  1. /*
  2. @Node Header
  3. @Next Macro
  4.  
  5. Code:       init.rexx
  6. Author:     Russell Leighton
  7. Revision:   16 Feb 1994
  8.  
  9. Comments:
  10.  
  11. This script sets up Ami2D on startup.  It is used by ami2d.rexx, but can be
  12. executed  to  reset to an initial state.  This script sets up Ami2D for use
  13. with structural analysis, plane strain problems.
  14.  
  15. @EndNode
  16. @Node Macro
  17. */
  18.  
  19. address ami2d
  20.  
  21. 'reset all'
  22.  
  23. 'set struct 0'
  24. 'set stress 1'
  25. 'set strain 2'
  26. 'set thermal 1'
  27.  
  28. /* iso(id,E,nu,alpha,Tref) */
  29.  
  30. 'set la=a2*a3/((1+a3)*(1-plane*a3))'
  31. 'set mu=a2/(2*(1+a3))'
  32. 'set et=a4*a2/(1-plane*a3)'
  33. 'alias iso=mat(a1,0,struct,la+2*mu,la+2*mu,la+2*mu,mu,la,la,0.0,la,0.0,0.0,thermal,et,et,a5)'
  34.  
  35. /* therm(id,k) */
  36.  
  37. 'alias therm=mat(a1,0,thermal,a2,0.0,0.0,a2,0.0,0.0,0.0,0.0,0.0,0.0)'
  38.  
  39. /* tri3(id,mid,n1,n2,n3) */
  40.  
  41. 'alias tri3=elem(a1,a2,3,1,a3,a4,a5)'
  42.  
  43. /* tri6(id,mid,n1,n2,n3,n4,n5,n6) */
  44.  
  45. 'alias tri6=elem(a1,a2,6,3,a3,a4,a5,a6,a7,a8)'
  46.  
  47. /* quad4(id,mid,n1,n2,n3,n4) */
  48.  
  49. 'alias quad4=elem(a1,a2,4,4,a3,a4,a5,a6)'
  50.  
  51. /* quad8(id,mid,n1,n2,n3,n4,n5,n6,n7,n8) */
  52.  
  53. 'alias quad8=elem(a1,a2,8,4,a3,a4,a5,a6,a7,a8,a9,a10)'
  54.  
  55. /* quad9(id,mid,n1,n2,n3,n4,n5,n6,n7,n8,n9) */
  56.  
  57. 'alias quad9=elem(a1,a2,9,9,a3,a4,a5,a6,a7,a8,a9,a10,a11)'
  58.  
  59. /* dispx(nid,value) */
  60.  
  61. 'alias dispx=pbc(a1,1,0,1,a2)'
  62.  
  63. /* dispy(nid,value) */
  64.  
  65. 'alias dispy=pbc(a1,1,0,2,a2)'
  66.  
  67. /* forcx(nid,value) */
  68.  
  69. 'alias forcx=pbc(a1,2,0,1,a2)'
  70.  
  71. /* forcy(nid,value) */
  72.  
  73. 'alias forcy=pbc(a1,2,0,2,a2)'
  74.  
  75. /* dispr(nid,value) */
  76.  
  77. 'alias dispr=pbc(a1,1,0,1,a2)'
  78.  
  79. /* dispz(nid,value) */
  80.  
  81. 'alias dispz=pbc(a1,1,0,2,a2)'
  82.  
  83. /* forcr(nid,value) */
  84.  
  85. 'alias forcr=pbc(a1,2,0,1,a2)'
  86.  
  87. /* forcz(nid,value) */
  88.  
  89. 'alias forcz=pbc(a1,2,0,2,a2)'
  90.  
  91. /* slope(eid,side,value) */
  92.  
  93. 'alias slope=dbc(a1,1,0,1,a2,a3)'
  94.  
  95. /* press(eid,side,value) */
  96.  
  97. 'alias press=dbc(a1,2,0,1,a2,a3)'
  98.  
  99. /* tract(eid,side,value) */
  100.  
  101. 'alias tract=dbc(a1,2,0,2,a2,a3)'
  102.  
  103. /* tempr(nid,value) */
  104.  
  105. 'alias tempr=pbc(a1,1,1,1,a2)'
  106.  
  107. 'set scale 1'
  108. 'set material 1'
  109. 'set disp 0'
  110. 'set force 0'
  111. 'set press 0'
  112. 'set tract 0'
  113.  
  114. 'set temp=u1'
  115.  
  116. call 'ami2d:rexx/pstrain'
  117.  
  118. exit
  119.